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Abstract 

"q^" ■ Some topics on heavy-ion collisions are reviewed with emphasis on those 

O . which are expected to be specially relevant at the Large Hadron Collider pro- 

^ ' gramme. 

1 Introduction 

r**- ! Heavy ions will be collided at the LHC with energies thirty times larger than ever before: each individual 

nucleon inside the colliding nuclei will be accelerated at three times the top Tevatron energy in the 
JZ: ■ centre-of-mass frame — which is already ten times more than that in AuAu collisions at the Relativistic 

Heavy Ion Collider (RHIC) in Brookhaven. The energy frontier reached in heavy-ion and proton-proton 
Oh' collisions is, for the first time, the same, since the same machine is used to accelerate the two different 
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systems 1 . Most of the currently targeted physics quests refer, however, to very different energy scales; 
those in heavy ions correspond to typical QCD scales (Aqcd < 1 GeV) while the scales for the principal 
LHC searches in the proton-proton programme lie 2-3 orders of magnitude above. 

The main goal of heavy-ion collisions is to form and characterize a macroscopic (in QCD scales) 
state of deconfined quarks and gluons in local thermal equilibrium. The large amount of energy liberated 
CN ■ in the collision is distributed over distances of the transverse size of an atomic nuclei {Ra ~ A 1 / 3 fm 

~ 6 fm for a lead nucleus) and expected to give rise to collective effects with distinctive signatures. 
This distribution of energy in the largest possible length scale is the most probable configuration in any 
■ hadronic collision. On the other hand, the configuration in which a large amount of energy is concentrated 

in a small region of phase space occurs very rarely, but interestingly leads to the largest probability of 
producing something unknown like a new particle. The second possibility is specially relevant for the 
searches in the LHC proton-proton programme. With the study of hard processes it is also becoming 
more and more important for the heavy-ion programme. 

In some cases, a nucleus-nucleus collision can be considered as a superposition of independent 
proton-proton collisions. This happens, for example, for hard enough probes, as photon production, in 
which the scales probed are much larger than those ruling the collective behaviour (or, in other words, 
those for which the production time is much shorter than the typical times involved in collective be- 
haviour). These probes are very important as benchmarks. In general, however, the output of a heavy- 
ion collision is very different from such a simple superposition. The study of processes sensitive to the 
degree of collectivity of the system is the main goal of heavy-ion collisions. For example, a powerful 
way of looking for collective behaviour is by studying the azimuthal distribution of particles produced at 
a given transverse momentum px- if the collisions were independent, the particles would be uniformly 
distributed (within statistical fluctuations) while the presence of azimuthal asymmetries would indicate 
that some non-trivial phenomenon is taking place. Different other probes of this collective behaviour 
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have been proposed both in the soft and the hard sector of QCD — ruled by non-perturbative and per- 
turbative dynamics, respectively. In the next sections we shall review some of these probes and the new 
opportunities opened at the LHC. 

In parallel to the PbPb programme, pPb collisions will also be made at the LHC. One of the main 
goals will be in this case to provide the essential benchmark for a correct interpretation of the hot matter 
probes — the energy densities reached in a pPb collision at the LHC are not expected to form a hot 
medium in the final state. For example, the knowledge of the parton distribution functions is deficient 
at small- ssj, those most relevant for the LHC, and especially for the case of gluons. Moreover, this 
type of collisions present an interesting physics case per se, where non-linear terms in the evolution 
equations are expected to become important, marking the beginning of a new regime of QCD where the 
parton distributions are saturated. Extra physics opportunities are the ultraperipheral collisions and the 
measurement of cross-sections of cosmic -ray interest. 



2 The LHC heavy-ion collider and the experiments 

At the LHC, the centre-of-mass energy of a collision of two different systems with charge Z\, Z2 and 
atomic numbers A%, A2 with Z = A = 1 for a proton is (see Ref. (H for an updated report on the 
subject of nuclear collisions at the LHC) 



Z1Z2 



^s NN ~ + , (1) 

where N N in the subindex refers to the energy per nucleon inside the colliding nucleus and ^fs^, is the 
corresponding energy in pp collisions. This gives the maximum energies 



\J s nn — 5.5 GeV for PbPb collisions, 
V s nn — 8.8 GeV for pPb collisions. 

Additionally, for non-symmetric systems — as in the projected pA collisions — there is a rapidity shift 

due to the fact that the centre-of-mass frame of the pA collision does not coincide with the laboratory 
centre-of-mass frame. This rapidity shift is Ay ~ 1 for proton-lead collisions and would need to be 
taken into account for the comparison with PbPb data. 

Another relevant quantity is the maximum luminosity for the heavy-ion programme which is esti- 
mated to be 

L NN ~ 10 27 omTV -1 (3) 

where NN refers again to luminosity per nucleon. This luminosity per nucleon is ~ 7 orders of mag- 
nitude smaller than the corresponding nominal luminosity for a typical pp run — L pp ~ 10 34 cm _2 s _1 . 
A typical year when the LHC reaches its design features will be eight months of pp run and one month 
dedicated to the heavy-ion programme (including nucleus-nucleus or proton-nucleus collisions at dif- 
ferent energies or pp at smaller energies for benchmarking in the PbPb runs). This reduces the integrated 
luminosity by another factor of ten with respect to the proton-proton programme. 

Heavy-ion collisions will be detected by ALICE [2l — the LHC experiment dedicate to that 
purpose — and both CMS [4] and ATLAS 0. The three experiments have complementary capabilities 
for the study of the different probes proposed. They are, for the same reason, flexible enough to look for 
unknown or not yet considered processes. 

Kinematical reaches are, at the LHC, orders of magnitude different than the present ones both both 
for large virtualities and small- x — see Fig. [Qfor the kinematical reach explored in a typical year of pPb 
collisions. With these numbers it becomes clear that the LHC will provide exploration of completely 
new regimes of QCD under extreme conditions. 
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ALICE+CMS reach in 1 yr. pA(Ap) collisions 




Fig. 1: (Left) Kinematical reach in one year of pPb collisions at the LHC. PbPb collisions will have a similar reach. 
(Right) Reach in transverse momentum (energy) for different processes in PbPb collisions at the LHC, figure from 
Ref. 0. 

3 QCD at high temperatures 

QCD is the theory of strong interactions: it describes the interactions between hadrons (protons, 7r's, 
etc.) which are the asymptotic states of the theory at small temperatures and densities. They form 
nuclear matter, which constitutes the majority of the observable mass of the Universe. Hadrons are 
colourless objects. The Lagrangian of QCD is, however, written in terms of quarks and gluons which are 
the building blocks of the hadrons — hadrons are composite particles. Quarks and gluons are colourful 
objects. Colour is the charge of QCD and the fact that the bosons carrying the interaction, the gluons, are 
also coloured allows them to interact. This makes the theory very different from its abelian counterpart, 
QED. Quarks are colour triplets corresponding to the fundamental representation of the QCD gauge 
theory group SU (3). There are six flavours of quarks with fractional electric charges and very different 
masses: 



charge = 2/3 


u (~ 5 MeV) 


c(~ 1.5 GeV) 


t (~ 175 GeV) 


charge = -1/3 


d (~ 10 MeV) 


s (~ 100 MeV) 


b (~ 5 GeV) 



One normally distinguishes between light (u, d, s) and heavy (c, b, t) quarks. Two extremely 
important properties of QCD are confinement and asymptotic freedom. The first corresponds to the non- 
existence of colourful asymptotic states — i.e., colour is confined in regions smaller than a typical QCD 
scale, Aq^ d ~ 1 fm. In particular, no free quarks or gluons can exist as asymptotic states. Asymptotic 
freedom, on the other hand, corresponds to the behaviour of the theory at small distances or large scales. 
The QCD coupling constant becomes smaller and smaller — the interaction disappears — at these large 
scales or small distances QUI. 



A simple picture for the confinement is that of a string with a quark and an antiquark at each end. 
The potential between the qq pair can be written as [9] (this is known as the Cornell potential): 

V( r ) ~ -— + ar. (4) 
r 

The linear part of the potential becomes relevant at large distances; a is called the string tension. So, 
if we try to pull apart the quark from the antiquark, the linear potential will make the amount of energy 
grow indefinitely; an isolated quark in the vacuum has infinite energy. This, however, does not actually 
happen since when the energy in the string is larger than m q + mq, the string can break by pair creation 
and form two different strings — two different hadrons. In the limit m g — ► oo the string cannot break. 

Another important property of QCD is chiral symmetry. In the absence of quark masses, the QCD 
Lagrangian splits into two independent quark sectors, 

£-QCD = £ g luons + iq~Ll^D^q L + iqR^D^q R . (5) 

This Lagrangian is symmetric under SU (N) l x SU(N)r with N the number of massless quarks. How- 
ever, this symmetry is not observed experimentally. The mechanism of spontaneous symmetry breaking 
does not need in QCD the introduction of a new field, at variance with the electro-weak sector, whose 
symmetry is broken by the presence of a Higgs field. In QCD, the vacuum |0) is not invariant under the 
transformation, the chiral condensate is non-zero 

(0\q L qR + q R qL\0) / 0. (6) 

The Goldstone theorem states that, when continuous symmetries are spontaneously broken, associated 
massless bosons should appear. In the case of N = 2 this corresponds in QCD to the three pions ir^, tt° 
whose masses are much smaller than other typical hadronic masses. The case for N = 3 includes the 
lowest-mass strange mesons. 

So, two main properties of the QCD vacuum are confinement and chiral symmetry breaking. A 
relevant question is then: Is there a regime where these symmetries are restored?. Intuitively, asymptotic 
freedom should lead to a deconfined free gas of quarks and gluons at very large temperatures and, indeed, 
this was proposed soon after asymptotic freedom was discovered iflOl 2 . The generic name for this new 
phase of matter is quark-gluon plasma (QGP). 

The asymptotic state of T — » oo is a gas of free quarks and gluons, but the situation at experimen- 
tally reachable temperatures is not so easy. The temperature at which the transition from a hadron gas 
to a quark-gluon plasma takes place, T c ~ 160 -j- 190 MeV 3 , is not large enough to allow the use of 
perturbative techniques at temperatures close to it. So, although the use of effective models is still a very 
active field of research, most of the information about the transition and the properties of the matter in its 
vicinity comes from lattice QCD calculations, see for example, Ref. lfl4l for a description of the QCD 
thermodynamics as studied by lattice. 

A first example of this collective behaviour is the equation of state. The pressure or the energy 
density measured in units of T 4 of an ideal gas is, according to the Stefan-Boltzmann law, proportional 
to the number of degrees of freedom in the system: so, for a free gas of pions this quantity is iVJ of = 3, 
while for a free gas of quarks and gluons this quantity is much larger, f = 2x2x3, N§ { = 2x8, 
counting spin, colour and (two) flavour states. This different behaviour is observed in lattice calculations 
where a jump at the transition temperature, T c , appears (Fig. |2]). Interestingly, the lattice results also 
signal a significant departure of the ideal gas behaviour, e = 3p, at temperatures close to T c (Fig. |2]). 

2 The original proposals of a different, deconfined, phase of hadronic matter at high temperatures is based also on the 
exponentially rising spectrum of hadronic states 1111 which predicted the existence of a limiting temperature, then identified as 
the transition temperature to the QGP. 

3 At present there is a theoretical debate on the actual value of this quantity as different lattice groups obtain results approxi- 
mately in the quoted range. The disagreement seems to be in how the continuum limit is taken, and the corresponding matching 
from lattice to physical units. See Ref. 1 12, 13] for details. 
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Fig. 2: Left: Entropy, s = e + p in units of T 3 versus temperature computed in lattice. Right: Trace anomaly in 
units of T 4 as a function of the temperature. Figures from Ref. ifTTl 



In order to understand the nature of the transition order parameters are needed. In QCD with 
massless quarks, the chiral condensate is the order parameter (see Ref. [ 15] for a recent study), while in 
the infinite mass limit the order parameter is the Polyakov Loop. The order of the transition depend on 
the actual value of the quark masses and most calculations agree in the absence of a real phase transition 
at zero baryochemical potential — the transition would just be a rapid cross-over. 

The behaviour of the potential between heavy quarks is also of interest for phenomenological 
applications. Although some discussion about the precise meaning and the definition of a potential exists 
in the case of a hot medium lfl6l . a screening leading to a non-confining potential is expected to appear 
at some point above the phase transition and, correspondingly, the heavy quark bound states would cease 
to exist if this temperature were reached. 

4 General properties of heavy-ion collisions 

Nuclei are very extended objects for all the scales of interest in high-energy physics. For this reason, 
the geometry of the collision plays a central role in the analysis and interpretation of the experimental 
results. In the centre-of-mass frame, owing to Lorentz contraction in the longitudinal direction, the two 
nuclei can be seen as two thin disks of transverse size 2Ra — 2A 1 / 3 fm. Some relevant quantities are 

1. The impact parameter which is the distance between the centres of the two colliding nuclei. The 
impact parameter characterizes the centrality of the collision. A central collision is one with small 
impact parameter in which the two nuclei collide almost head-on; a peripheral collision, on the 
contrary, is one with large impact parameter. A special case is that of ultraperipheral collisions in 
which the two nuclei interact only through the photons created by the large electromagnetic field 
of the ions lfT8l . 

2. The number of participant nucleons, A^ part , within the colliding nuclei which is the total number 
of protons and neutrons which take part in the collision. The rest are sometimes called spectators 
and continue travelling almost unaffected; they are normally measured in forward detectors such 
as the Zero-Degree Calorimeter (ZDC) to help in fixing the centrality of the collision. 

3. The number of collisions, iV co ii, is the total number of incoherent nucleon-nucleon collisions. 



4.1 The Glauber model and the geometry of the collision 

The usual way of computing the geometrical quantities enumerated above is by a probabilistic model 
due to Glauber QjQ. The starting point is the density of nucleons in a nucleus A, p(z, b), for a given 
longitudinal z and transverse b positions. The nuclear profile function or thickness function is the only 



relevant quantity in most of the calculations and is denned as 

poo 

T A (b) = / dzp(z,b). 



(V) 



With the normalization f dbT^b) = 1, the individual probability of a nucleon-nucleon interaction at 
a given impact parameter is T A (b)a^ff. If we now consider a proton-nucleus collision, the probability 
that the proton interacts with n nucleons inside the nucleus is simply given by 



P(n,b) 



A 



n 



i-T A {h)^i 



A—n 



T A (b)a^% 



(8) 



from which we can compute the number of collisions for a given impact parameter as 

A 



<n(b) = < ar t(b) = J>P(n,b) = AT A (b)a%% , 



(9) 



n=0 



where the number of participants of the nucleus A, N^ Art , coincides in this case with the number of 
collisions N co \\ — each nucleon in nucleus A interacts only once. The total number of participants 
including the colliding proton is then + 1. The inelastic cross-section can be computed as 



<A = / db£i>(n,b) 



n=l 



db 



1 - T A {b)o%$ 











~ j db 


1 — exp 



-AT A (b)a%% 



(10) 

where the limit of large A has been taken for the last equality. For a nucleus-nucleus collision, the 
individual probability of nucleon-nucleon interaction at impact parameter b is within the optical ap- 
proximation — see for example I1201I2T1 for experimental applications. 



dsT A (b)T B (b - s)a i N e N = T AB (b)a%% . 



inel 



(11) 



T A b (b) is also known as the nuclear overlap function. The corresponding probability of n NN collisions 
is 



P(n,b) 



AB 

n 



1 - 2AB(b)oAW 



AB—n 



T AB (b)a 



which gives the inelastic cross-section and number of collisions as 



<f(b) 



db 



1- [l-T AB (b)a^ 











^ AB' 


~ J 


j db 


1 — exp 



inel 
NN 



-ABT AB (b)a%i 



(12) 

,(13) 
(14) 



n=0 



The number of participants of nucleus A for a given impact parameter is given by the generalization of 
(©; for this we have to single out one nucleon in nucleus B at transverse position b by writing BT B (b) 



< rt (b) 



dsBT B (s)aZ\b-s) = dsBT B (s)exp —AT A (b — s)a 



'NN 



(15) 



With a similar expression for the number of participants in nucleus B, to give the total number of partic- 
ipants as 

AW(b) = JV£rt(b) + iV p B art (b) . (16) 
In Fig. [3] a picture of the geometry described above can be found. The measure of the geometry in 
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Fig. 3: Picture for the geometry of a high-energy heavy-ion collision. 

heavy-ion collisions is performed using several methods, with he possibility to cross-check. For example, 
knowing the number of particles produced in each elementary NN collision, the multiplicity distributions 
can be related to the number of participants or collisions. Another method is by measuring with ZDC 
the number of neutrons which did not interact (spectators). In this case, the number of participants is 
estimated by the total number of nucleons minus the result from the ZDC. As the impact parameter is not 
a quantity of direct experimental access, a usual way of quoting the centrality of a class of measurements 
is by a percentage, e.g. 0-5 % refers to the 5% most central collisions; 0-50% to the 50% most central 
collisions; 10-20% to the 20% most central collisions excluding the 10% most central ones and so on. 

4.2 The multiplicities 

The total number of particles measured in the detectors after a heavy-ion collision (multiplicity) gives 
information about the energy density reached. It also gives information about the centrality and, in 
general, about the main global properties of the created medium. The knowledge of this quantity and its 
corresponding theoretical description is, then, one main issue in heavy-ion collisions. 

In hadronic collisions, most of the particles are produced at small transverse momentum, belonging 
to the soft sector of QCD in which the use of perturbative techniques is not possible. This makes the 
theoretical description difficult and in most of the cases models are employed. Specially interesting is 
the formalism based on saturation of partonic densities in the initial nucleus wave function which allows 
weak coupling calculations to be made in a regime of large occupation numbers. This formalism is 
generically known as the Colour Glass Condensate. 

The most naive estimate of the multiplicity in AA collisions would be that of the corresponding 
proton-proton collision at the same energy and multiplied by the number of incoherent AW collisions. 
This number, however, already grossly overestimates the number of produced particles at RHIC, demon- 
strating the presence of collective behaviour. All the models of multiparticle production prior to RHIC 
data ll22l that did not include any type of collectivity failed to reproduce the data. The surviving models 
present, however, a quite large uncertainty in the extrapolation to the LHC, see Fig. |4] In the rest of the 
section we present some results based on saturation physics. The interested reader can find the relevant 
references to the different models in the workshop |2"4l . 
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Charged multiplicity for r|=0 in central Pb+Pb at ^1^=5.5 TeV 
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Fig. 4: Different predictions for the charged particle multiplicities at central rapidities at the LHC. Figure from 
Ref.j23); results are a compilation from Ref. [24] . 



4.2.1 Ideas based on saturation ofpartonic densities 

The density of gluons 4 grows very fast at small x. The intuitive picture of saturation of partonic densities 
is that of overlapping partons, mainly gluons, in transverse plane due to this growth (Fig. [5]>. In the 
infinite momentum frame, in which the hadron/nuclei is moving asymptotically fast, the colliding object 
is a Lorentz-contracted disk of radius Ra ~ A 1 / 3 . If we do DIS with this object, a virtual photon with 
virtuality Q S> Aqcd will scatter incoherently with the partons inside it: counting the number of partons 
with transverse size ~ l/Q. Owing to partem multiplication the number of gluons in the proton/nucleus 
grows very fast with decreasing x, i.e., with increasing energy. At some point, the transverse size occu- 
pied by the partons is of the same order as the total transverse size of the hadron. This overlap makes 
gluon fusion probable, with probability 0(a s ), and the picture of incoherent scatterings breaks down. 
This defines the saturation scale which can be geometrically estimated as 

a s -^Axg(x,Q 2 SSLt ) ~ ttR 2 a . (17) 

^sat 

For a gluon density which grows as xg(x, Q 2 ) ~ x~ x , the behaviour of the saturation scale as given by 
this naive geometrical estimate is 

Q 2 sat ~ x-W . (18) 

Here, the gluon distribution of bound nucleons is taken as Axg(x, Q 2 ) — there is no shadowing see 
below — and the nuclear radius Ra — A 1 ' 3 . 

Intuitively, when Q SS/t 3> Aqcd this should be the relevant scale in the problem and the parton 
distributions should depend only on this scale. Indeed, one of the predictions from the Colour Glass 

4 The number of gluons is not a well defined quantity in a field theory like QCD. All the definitions in this section can be 
done in a more strict way, e.g. by using scattering amplitudes, but this is a useful intuitive picture which contains most of the 
physics involved. 
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Fig. 5: Saturation of partonic densities: an intuitive picture of saturated gluon distributions in the transverse plane. 
Each gluon is depicted by a small disk whose number grows with atomic number and/or energy of the collision. 



Condensate is the scaling of the parton distribution functions with the saturation scale Q 2 at (geometric 
scaling). This scale encodes the geometry of the colliding object as well as the dynamics of the gluon 
multiplication and fusion. 

These ideas lead to a very developed formalism aimed at describing, not only the hadron/nucleus 
wave function, but also other aspects of QCD at asymptotic energies like multiparticle production or 
thermalization in heavy-ion collisions, see, for example Refs. I25ll26ll27ll28l [291 for lectures on the topic. 
Here we present for illustration a simple implementation of this physics based solely on experimental 
data and the assumption of geometric scaling ll30l . 

i) Geometric scaling in lepton— proton and lepton— nucleus collisions 

All data for the photoabsorption cross-section <r 7 * p (x, Q 2 ) in lepton-proton scattering with x < 
0.01 have been found [31] to lie on a single curve when plotted against the variable Q 2 /Q 2 p , with 
Q 2 ~ x~ x and A ~ 0.3. Geometric scaling is usually motivated in the QCD dipole model ll32l where 
the total j*h cross-section reads 

4 h L {x,Q 2 ) = J drj\z\^ L {Q 2 ,r,z)\ 2 a^ p {r,x). (19) 

Here ^t,l are the perturbatively computed transverse and longitudinal wave functions for the splitting 
of 7* into a qq dipole of transverse size r with light-cone fractions z and (1 — z) carried by the quark and 
antiquark respectively. Both for a proton [h = p] and for a nucleus [h = A], o"j ip (r, x) can be written as 
an integral of the dipole scattering amplitude Nh over the impact parameter b, 



a£ p (r, x) = 2 J db N h (r, x; b) . (20) 

In this setting, geometric scaling corresponds to the condition iV^(r, x;h) = Nh(rQ s ^(x,b)). This 
can be seen by rescaling the impact parameter in (1201 in terms of the radius of the hadronic target, 



b = h/^R 2 , 

a^(x,Q 2 )=7TR 2 J dr^dz\^ L (Q 2 ,r,z)\ 2 2 J db N h (rQ s , h (x, b)) . (21) 

For a trivial impact-parameter dependence of the saturation scale, Q Sj h(x,b) = Q St h{x) @(Rh — b), 
and since L (Q 2 , r, z)\ 2 is proportional to Q 2 times a function of r 2 Q 2 , Eq. (|2TI ) depends solely 
on Th = Q 2 /Q 2 h (x). For realistic functional shapes of the form Q S: h(x,b) oc f(b/Rh), the same 
r-dependence results if Q 2 h (x) is defined as an appropriate b-average of Q 2 h (x,b). In the case of 

7* A interactions, geometric scaling is the property that the yl-dependence of the ratio Oj, L /ttR\ can be 
absorbed in the A-dependence of this impact-parameter independent saturation scale Q Sj \(x), 

^* A (t A ) _ qf*{T A ) _ 
kR 2 a ttR 2 



Now we need to fix the yl-dependence of the saturation scale, for which we make the assumption — see 



where the nuclear radius is given by the usual parametrization Ra = (1.12A 1 / 3 — 0.86J4" 1 / 3 ) fm. We 
treat 5 and ttR 2 as free parameters to be fixed by data. 

In Fig. [6] we plot the experimental ^y*p data [33 ] with x < 0.01 as a function of r p = Q 2 /Q 2 p . 
For Qs p , we use in this plot the Golec-Biernat and Wiisthoff (GBW) parametrization P4l with Q 2 p = 
(x/x )- x in GeV 2 , x = 3.04 • 10" 4 and A = 0.288. In order to proceed, we need a parametrization 
of the experimental data (a functional form for the scaling curve). We find that the data 11331 are well 
parametrized by 

a^P(x,Q 2 ) = $(r p ) =<7 [7£ + r(0,0+ln£] , (24) 

where is the Euler constant, T (0, £) the incomplete F function, and £ = a/Vf, with a = 1.868 and 
b = 0.746. The normalization is fixed by oq = 40.56 mb. 

To determine Q 2 A , we compare the functional shape of ((24l to the available experimental data 
for j*A collisions with x < 0.0175 051 l36l l37l. using £ = a/r b A . The parameters 5 and ttR 2 in 
Eqs. (I22l-(l24l are fitted by \ 2 minimization adding the statistical and systematic errors in quadrature. 
The data sets [35J, Il36l . and ll37l have additional normalization errors of 0.4%, 0.2%, and 0.15%; the 
quality of the fit improves by multiplying the data by the factors 1.004, 1.002, and 0.9985, respectively. 
We obtain 5 = 0.79 ± 0.02 and ttR 2 , = 1.55 ± 0.02 fm 2 for a x 2 /dof = 0.95 — see Fig. [6] for 
comparison. If the normalizations are all set to 1, we obtain an almost identical fit with 5 = 0.80 ± 0.02 
and ttR 2 = 1.57 ± 0.02 fm 2 for a x 2 /dof = 1.02. If we impose 5 = 1 in the fit, which corresponds 
to Q 2 A oc A 1 / 3 for large nuclei, a much worse value of x 2 /dof = 2.35 is obtained. We conclude 
that the small-x experimental data on ^y*A collisions favor an increase of Q 2 sA faster than A 1 / 3 . The 
numerical coincidence b ~ 5 is consistent with the absence of shadowing in nuclear parton distributions 
at Q 2 > Q 2 s A . 

ii) Geometric scaling and multiplicities in heavy-ion collisions 

The fact that geometric scaling is found both in lepton-proton and -nucleus data is interesting, but 
we can extend the formulation a bit further and consider the case of multiplicities in symmetric hadronic 
or heavy-ion collisions — this means the same colliding systems and measure at central rapidities. This 
allows us to make several simplifications as only one saturation scale is involved. For the discussion, 
consider for simplicity a model of multiparticle production in which a factorization between the gluon 
distributions <fi(y, k, b) of the two nuclei A and B appears ll38ll 

dN AB 



dydptdh 



oc / dk <t> A (y, k 2 , b) B (y, (k - Pi ) 2 , b) , (25) 
Pt J 



where 0^(y, k, b) = / dr expjzr • k} Nh(r, x; b)/(27rr 2 ) l39l . y = Inl/x. For geometric scaling, 
<Pa(u, k 2 , b)= 4>(k 2 /Q 2 A (y, b)), we find the dependence 



dN AA 



dy 



y~0 



p? V QIa \ Q 2 



s.A 



^drdh^^dr-s) 2 ). (26) 
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Fig. 6: Geometric scaling for j*p (upper panel, data from Ref. If33l ). (middle panel, data from Refs. Il35l[36l0 
and the ratio of data for 7* A over the prediction from d24l i (lower panel). As an additional check, the lower plot 
also shows data for normalized with respect to 7*C ll37l . and divided by the corresponding prediction from 
Eq. (El. 



dy 



The integral is then independent of the colliding objects — geometric scaling has been used — resulting 
in 

ocQ s 2 A vr^. (27) 

Similar expressions arise in different models of hadroproduction If38l [39l l40l [4TTl . It is worth noting that 
factorization is not actually needed in (|25T ). any integrand with (£;/Q Sj A)-scaling leads to Eq. (|27"T ). see 
Refs. BT1I401 . In all these models, the hadron yield is assumed to be proportional to the yield of produced 
partons. 

Equation (l27"T ) relates the energy and nuclear size dependence of the multiplicities to the corre- 
sponding quantities measured in lepton-hadron collisions. In particular, the energy dependence is given 
by the GBW parameter A = 0.288. For the centrality dependence of (|27"1) , we use the known propor- 
tionality in symmetric A+A collisions between the number iV par t of participant nucleons and the nuclear 
size A. With Q 2 S A oc A 1 / 35 , and 5 = 0.79 ± 0.02, we thus obtain 5 , 



2 dN AA 



iVpart drj 



-A l ~ s 



NoV^ Xlrt • (28) 



part 

»7~0 



As seen in Fig. 13 this ansatz accounts for experimental data from the PHOBOS Collaboration ll42l on 
charged multiplicities in Au+Au collisions at yfs = 19.6, 62, 130, and 200 GeV/A. Even the p+p data 
( 11431 , as quoted in Ref. Ell) at ^fs = 19.6 and 200 GeV are accounted for by Eq. d28l) . Since all data are 
at mid-rapidity, the Jacobian between rapidity y and pseudorapidity rj is approximately constant. It has 



3 Notice a factor of 2 missing in the original paper [ 30 1. 
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Fig. 7: Energy and centrality dependence of the multiplicity of charged particles in Au+Au collisions (f28b com- 
pared to PHOBOS data |42|. Also shown in the lower panel are the p+p data |j43l and results for y/s — 5500 
GeV/A. 

been absorbed in the overall normalization Nq = 0.47 which is independent of the energy and the cen- 
trality of the collision. Figure |7] also shows the prediction of (l28T ) for LHC energy (^/s = 5500 GeV/A) 
— this result is quoted in Fig. @]as "Armesto et al" — and for smaller colliding nuclei. Equation (|28T ) 
implies that the energy and the centrality dependence of the multiplicity factorize, in agreement with the 
results by PHOBOS 1421 . 

Hi) Multiplicities computed from the QCD evolution equations 

Although the CGC approach is a very advanced theoretical construction, its implementation in 
actual phenomenological analysis is difficult, due, in part to the complexity of the equations to be solved 
and in part to unknowns in, for example, the initial conditions for evolution equations, etc. For this 
reason, most of the analyses of the data in the last years made strong simplifications or, simply, tried to 
find generic features in the data expected from the CGC such as the geometric scaling described above. 
Another difficulty was the use of leading order results which are known to produce too fast evolution 
of the parton distributions. A breakthrough in this direction is the computation fl4l l45l l46l of part of 
the non-leading corrections to the equations, leading to dependences in agreement with expectations 
from data. These equations are solved numerically in Ref. [47] to describe the energy and rapidity 
dependences of the multiplicities in nuclear collisions for the first time using first-principle calculations 
of the evolution equations. The results of this study are plotted in Fig. [8] and quoted in Fig. @] as 
"Albacete". 



5 Soft probes and the bulk 

The soft part of the produced spectrum provide very useful tools to characterize the collective properties 
of the nuclear collisions — and eventually the properties of a new state of matter. It is not our purpose 
here to give an extensive review of all probes proposed in heavy-ion collisions. We will, instead, limit 
ourselves to the ones specially relevant in the last years, with the advent of collider energies in nuclear 
collisions at RHIC. 

A powerful way of looking at the degree of collectivity reached in a heavy-ion collision is by 
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Fig. 8: Rapidity and energy dependence of the multiplicity of charged particles in Au+Au collisions (|28T > compared 
to PHOBOS data |42| and the corresponding extrapolation to LHC energies. 



studying the azimuthal distribution of particles produced at a given transverse momentum p?' if the col- 
lisions were all independent, the particles would be uniformly distributed (within statistical fluctuations) 
while the presence of azimuthal anisotropics would indicate that non-trivial phenomena are taking place. 
The hydrodynamic modeling of the heavy-ion collisions attempts to describe such phenomena. 

5.1 The hydrodynamical evolution 

The use of hydrodynamics to solve the space-time evolution of hadronic collisions goes back to Landau 
in the 1950s l49l and has been established as one of the main fields of research in heavy-ion collisions. 
In those dynamical situations in which the mean free path of the particles in the medium is very small, 
medium properties such as energy density, pressure, or temperature can be described in a hydrodynamical 
approach. The basic equation corresponds to the conservation of the energy-momentum tensor 

= _ (29) 

Neglecting viscosity, the energy-momentum tensor can be written as 

T^ v = (e + p)u»u u - pg» u , (30) 

where e, p, and are the (local) energy density, pressure, and four-velocity of a fluid element. The 
latter is normalized as u^u^ = 1. Lectures on hydrodynamical modelling of heavy-ion collisions can be 
found, for example, in Refs. ll50"ll5Tll . 

5.1.1 The Bjorken model 

A very simple model for the hydrodynamical description of a heavy-ion collision due to Bjorken Il52l 
captures some of the properties of the evolution, neglecting all the transverse dynamics: Assume the col- 
lision of two big nuclei at very high energy so that, owing to Lorentz contraction, they can be considered 
as two disks of infinite transverse size. Define the spatial rapidity and the proper times as 

y= l \ g t -^ ; T=Vt 2 + Z*. (31) 

2 t — z 
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At asymptotic collision energies, boost invariance is a good approximation if we are interested in the 
central rapidity region. To see this, recall that a boost corresponds to an additive term in the rapidity, 
y' = y + yboost ; if the energy of the colliding objects is very large, small longitudinal boosts (compared 
with the collision energy) should not affect the results, i.e., the result does not depend on the rapidity. In 
this case, the initial conditions are 

p(r); e(r); = 7 (1, 0, 0, z/t) , (32) 

and the hydrodynamic equations simplify to 

^ + ^ = 0. (33) 
dr t 

We need an equation of state, relating e and p, to close the equations. For an ideal EoS, e = 3p, the 
solutions are 

< T ) = ^> T (r) = ^k, (34) 

which tells how the energy density or the temperature decrease with proper time due to longitudinal flow 
in a boost-invariant ideal fluid. 



5.2 Transverse flow produced by hydrodynamics 

The above model is a very simplified version of the hydrodynamical evolution and, in particular, one very 
important feature of realistic heavy-ion collisions, the presence of transverse expansion, is not described. 
The sophisticated hydrodynamical calculations used nowadays include all these features, but for our 
purpose it is instructive to consider the Euler equation, contained in Eq. 



f— fv,. <3 5 ) 
dt e + p 

Equation (T331) tells us that gradients of pressure produce acceleration of fluid elements. This is especially 
relevant for the transverse plane, where, taking again the ideal case p = e/3, gradients of energy density 
are defined by the geometry of the collision. These transverse energy density profiles translate into 
a transverse expansion of the medium. The most interesting case is that of non-central collisions in 
which the medium has a characteristic almond shape with a strong azimuthal asymmetry in the energy 
density profile. This spatial asymmetry leads, hence, to an azimuthal-dependent transverse flow, i.e., the 
momentum of the particles after the freeze-out is strongly dependent on the azimuthal angle. In Fig. [9] 
an intuitive picture of the effect is depicted: i) two nuclei collide semiperipherally and the overlapping 
area is non-symmetric in the transverse plane; ii) the gradient of energy density is larger in the direction 
labelled as x than in the one labelled as y — these directions are normally called in-plane and out- 
of -plane for a plane, reaction plane, defined by the two beam directions; iii) this asymmetry leads to 
different density and pressure fields, producing azimuthally dependent accelerations, maximal in the y 
direction and minimal in the x direction. 

The transverse flow is one of the most important measurements in favor of the hydrodynamical 
behaviour in heavy-ion collisions and it is normally parametrized by the first non-trivial coefficient (at 
mid-rapidity) of the Fourier expansion in the azimuthal angle, called V2, 

dN 

— oc l + 2v 2 cos(2(/>). (36) 
d(p 

A value of v 2 7^ signals the presence of azimuthal asymmetries in the final state. Note that for a 
finite number of particles (N) these asymmetries are always present, in an event-by-event study, due 
to limited statistics; in the extreme case of two particles in the final state the asymmetry is maximized 
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Fig. 9: Hydrodynamics of a non-central nucleus-nucleus collision in the transverse plane. 
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Fig. 10: Comparison of a hydrodynamical calculation ll53l with experimental data from [48]. The effects of 
different values for the viscosity and the initial conditions in the hydrodynamical evolution are studied. 



as energy-momentum conservation implies they are back-to-back. These statistical asymmetries are, 
however, reduced as 1/ y/N. 

At present the use of hydrodynamical models is one of the most active subject areas in heavy-ion 
collisions. State-of-the-art calculations are beginning to include viscous corrections to the numerical 
simulations. In Fig. [10] we present the comparison with data of one from these studies |[53ll . 

6 Hard Probes 

The study of processes with large virtualities involved, such as production of particles with large masses 
and/or large transverse momentum, has become one of the main fields of research in heavy-ion collisions. 
In these hard processes different scales appear and can be separated into the different contributions to the 
factorized cross-section, which can be symbolically written as 

a AB^h = f A ( Xli c?) ® f B (x 2 , Q 2 ) ® a(xi,x 2 , Q 2 ) ® A-h(*, Q 2 ) ■ (37) 



IS 



Here, the large virtuality, Q 2 , allows us to compute the partonic cross-section, a(xi,x 2 , Q 2 ), as a per- 
turbative expansion in powers of a s (Q 2 ). This perturbative partonic process takes place in a short time 
~ 1/Q, i.e., it describes the processes associated to short-distances. In collisions involving hadrons, there 
is always a non-perturbative contribution associated to long-distances (small scales O(Aqcd)) such as 
the size of the hadron. Two of them appear in (I37T ): on the one hand, the partem distribution functions 
(PDF) /a {x, Q 2 ), encoding the partonic structure of the colliding objects at a given fraction of momen- 
tum x and virtuality Q; and the fragmentation functions (FF), D(z,Q 2 ), describing the hadronization 
of the partem i into a final hadron h with a fraction of momentum z. In the nuclear case, these are the 
quantities which are modified when the extension of the colliding system interferes with the dynamics, 
while the short-distance part is expected to remain unchanged if the virtuality is large enough. These 
modifications could involve the non-perturbative initial condition as well as the evolution equations. In 
the latter case, non-linear terms become important. 

At the LHC, where proton-proton, proton-nucleus and nucleus-nucleus collisions will occur, a 
good knowledge of the PDFs is essential: proton PDFs are needed to compute the cross-sections of 
processes which are background for searches such as the Higgs, SUSY, etc.; nuclear PDFs are also 
needed for benchmarking those effects in nucleus-nucleus collisions not related with the production of 
the hot medium under study. 

Although the presence of the non-perturbative PDFs is unavoidable at the LHC, there are processes 
in which the long-distance term after the hard collision, the FF in Eq. (|37T >. is not present. Some examples 
are photon, Z, or production which are produced directly in the hard scattering. These processes, 
well described by perturbative techniques, also have an important role as benchmarks in heavy-ion colli- 
sions. 

The employement of hard processes as tools to characterize the medium properties uses the fact 
that the non-perturbative, long-distance, terms in (I37T ) are affected by the presence of a macroscopic 
medium. A conceptually simple example is the case of the J/^f, whose production cross-section can be 
written as 

* hh -* J/ * = fi(x h Q 2 )®f j (x 2 ,Q 2 )®a i ^(x u x 2 ,Q 2 )(0([ce] - J/*)) , (38) 

where now (0([cc] — > J /^)) is a purely non-perturbative quantity describing the hadronization of a cc 
pair in a given state (for example a colour octet) into a final J/*. In the case where a hot medium is 
formed in the collision we can expect that the potential between the heavy-quarks will be screened and 
hence the bound states cannot survive 11541 . Ideally this would correspond to a strong reduction of the 
hadronization term (C([cc] — > J leading to a suppression of the yield of J/^'s if the hot medium 
is created. In practice, however, the theoretical control over the different terms in (|49l ) and their nuclear 
modifications is not good enough and the interpretation of the experimental results is still not clear. 

Another well-known experimental example of these modifications is the strong suppression of 
particles produced at large transverse momentum in heavy-ion collisions Il55ll56l . The effects due to the 
presence of a hot medium in the structure of the high-p^ part of the produced spectrum is known under 
the generic name of jet quenching, see, for example, Refs. Il57ll58ll59ll60l . 

In the next sections we report on the present knowledge of some of these hard processes see, for 
example, ED |62j |63] |64) for additional experimental results. 

6.1 Nuclear parton distribution functions 

It is a well-stablished fact that the partonic structure of high-energy nuclei is different from the incoherent 
superposition of the component nucleons, see Ref. [65 ] for a recent review. Indeed, experiments of DIS 
with nuclei have measured different nuclear effects in the whole range of x which can be associated to 
modifications of the nuclear PDFs when compared with the free proton ones. In particular, it has been 
customary to give different names to these modifications depending on the relevant range of x under 
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consideration: 1) shadowing for the suppression observed at small (x < 0.05); 2) antishadowing for the 
enhancement at moderate values of 0.05 < x < 0.3; 3) EMC effect for the suppression observed in the 
region 0.3 < x < 0.7; and 4) Fermi motion for the enhancement when x — > 1. 

The observed effects on the structure functions measured in DIS with nuclei can be implemented 
as modifications of the non-perturbative input in a DGLAP analysis of the nuclear PDFs when compared 
with the free proton case. This type of analyses has reached a degree of sophistication similar to the ones 
for free protons, with studies done at LO and NLO and error analysis following the Hessian method. 
However, the small amount of experimental data, especially the region of small-x, makes the knowledge 
of the nPDFs not as precise as will be necessary for a correct interpretation of all the experimental data 
at the LHC. This deficiency will be cured at the LHC with the presence of a parallel proton-nucleus 
programme. 

6.1.1 Global fits 

As in the case of the free proton, the sets of nuclear PDFs are obtained by global fits to different sets of 
experimental data (Sll67llMl|63|73inil72l|73l|7l|7a. Nowadays, all groups include data on DIS with 
nuclei and Drell-Yan production in proton-nucleus collisions. Data on inclusive particle production in 
deuterium-gold collisions measured at RHIC was first included in the analysis of Ref. 11741 . showing, in 
particular, the feasibility of these analyses and the uncertainties to be faced when pA data from the LHC 
becomes available. 

Most of the groups perform the analysis at the level of the ratios 

defined as the ratio of the PDF for a proton inside the nucleus, f^(x, Q 2 ), to the corresponding one for 
a free proton, fi(x,Q 2 ). Assuming a given functional form for the ratios 091 ) at an initial scale Qq, their 
parameters are obtained by evolving both the numerator and the denominator by the DGLAP evolution 
equations [76] and fitting the used sets of data. 

The procedure to perform nPDF global fits follows closely that of the free proton case: 

1. Take a functional form for the ratios at the initial scale Qq ~ 1 GeV 2 with a number of free 
parameters {z}: Ri(A, x, Qq, {z}), i = g, u, u, d, d, s . . . 

2. Assuming a known set of free proton PDFs, evolve the numerator and denominator of d39l) sepa- 
rately to a larger scale Q 2 using the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equa- 
tions 11761 . 

3. Compute the theoretical values corresponding to each data set. This gives a value of x 2 for a 
particular choice of the parameters {z} 

N 

4. If the function defined by the parameters {z} is the minimum of \ 2 this is the final answer, other- 
wise, choose another {z} and start again from point 2. 

The main difficulty in the global fit analyses is to find an initial condition for the evolution (point 1. 
above) which is flexible to avoid biasing the result as much as possible but with a small enough number 
of free parameters so that a minimization procedure can be performed. This problem is common to the 
free proton analyses, but more severe for nuclei due to the much more limited number of experimental 
data available, specially in the small-x region. In practice, some assumptions are imposed for those 
kinematical regions and/or flavours less constrained by the data. For example, a common assumption is to 
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Fig. 11: Comparison of the average valence and sea quark, and gluon modifications at Q 2 = 1.69 GeV 2 and 
Q 2 = 100 GeV 2 for Pb nucleus from the NLO global DGLAP analyses HKN07 d, nDS CD, and EPS09NLO 
ll75l including the error estimates. Figure from Ref. f75l . 



take for the valence quarks R^ y (x, Qq) = R§ (x, Qq) = Ry(x, Qq) or for the sea quarks R^(x, Qq) = 
R^(x,Qq) = Rf(x,Qo) = Rs(x,Qq) 6 . In general, the different analyses give similar results for 
the ratios Rf(x, Q 2 ) for those regions in which experimental data exist, but differ otherwise. In Fig. 
[TT1 different sets for the ratios of PDFs of protons inside a lead nucleus over those in a free proton are 
plotted. The ratios of the different groups are rather similar for valence quarks: well constrained by 
DIS data except at small-x; and sea quarks: constrained by DIS and DY data except for very small and 
large-x. The case of the gluons is, however, less satisfactory as very different parametrizations (specially 
in the small- x region) lead to similar descriptions of the available DIS and DY data. In order to improve 
the situation, additional experimental data able to constrain the gluon distributions need to be included in 
the analyses. For the first time, dAu data on inclusive high-px particle production from RHIC have been 
included in Ref. f74ll . In particular, the strong shadowing of gluons found in this analysis can be traced 
back to the strong suppression of the forward rapidity negative particle yields measured by BRAHMS. 
It is worth emphasizing, however, that although the ratios are very different for the lower scales, they 
become more and more similar with evolution. In particular, the nuclear effects for the gluons at a 
virtuality of ~ 100 GeV are very similar for all the groups as showed in Fig. [12] 

These global fit analyses allows one also to compute the uncertainties for a given observable 
associated the nPDFs. The procedure, adapted from the free proton analyses, consists in computing the 
eigenvalues corresponding to the Hessian matrix H defined as 



TT 1 

l,J 2 dz,d Zi 



(41) 

_ o. ._ o 

a - i ' j j 



"3 

where {z®} corresponds to the values of the parameters which minimize x 2 - For N free parameters, the 
same number of eigenvalues {z{\ is obtained, each one defining one eigenvector. These eigenvectors are 



6 These assumptions are only taken at the initial scale Qq, but the equalities do not hold in general for larger virtualities as 
the evolution of each flavour is different. 
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computed by setting all the parameters {z{\ in the new basis to the values at the minimum {z®} except 
one, z~k- In order to compute the uncertainty, two new sets of nPDFs are then computed for = z9 + dz^ 
and Zk = z® — Sz^, where 5zk depends on the confidence level required. In this way, 2N different sets 
of nPDFs are defined. To compute the error band associated to a given observable, for example the jet 
cross-section in heavy-ion collisions at the LHC, one would compute this observable for the 2N different 
sets of nPDFs and the results added in quadrature. A computer code for practical applications including 
the nPDF error estimates is available in Ref. 1751 . 
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Fig. 12: Evolution of the gluon ratios Rf(x, Q 2 ) for the parametrizations of the different groups, EKS98, EPS08, 
HKN, nDS. The differences, present at the smallest virtualities, quickly disappear with increasing Q 2 . 



6.2 Non-linear evolution: saturation of partonic densities 

The saturation of partonic densities described in Section 14.2.11 finds its natural implementation in the 
presence of non-linear terms in the evolution equations. The first non-linear terms computed were the 
leading higher-twist corrections to the DGLAP evolution equations |[38l 1771 



dxg(x,Q 2 ) dxg(x,Q 2 ) 



9tt a 2 f 1 dy 1 
DGLAP ^ Q 



d log Q 2 d log Q 2 

where the first term is the standard DGLAP result f76l . linear in the PDFs, and R ~ 1 fm, the radius 
of the proton. The approach to the problem is, however, different nowadays and a generalization of 
the BFKL equation f79l 1801 is normally used in phenomenological approaches 7 instead of the DGLAP 
equation. Several reasons exist for this preference, in particular, the small- x evolution is theoretically 
enhanced by terms in a s log x, which are the ones resummed by BFKL. 

The most widely used evolution equation in phenomenological analyses is the Balitsky-Kovchegov 
(BK) ll8Tll82l equation. It gives the evolution with rapidity Y = In (s/so) = m (xq/x) of the scattering 
probability N(x, y, Y) of a qq dipole with a hadronic target, where x (y) is the position of the q (q) in 
transverse space with respect to the centre of the target. Defining 

f=x-y, 6 = x + |, fx = x — z, f 2 = y-z, (43) 

and neglecting the impact parameter dependence (which is justified for r <C b), the BK equation reads 

(r = \f\) 

dN(r,Y) a s N c f d 2 z r 2 



dY 



N f d z z r l 

— hr — i N ( r ^ Y ) + N ^ Y ) ~ N ^ Y ) ~ N(rx,Y)N(r 2 , Y)] . (44) 

7T J ZTT 



The coupling constant is fixed and the kernel is conformally invariant — this means, in particular that 
the solutions present geometric scaling. The NLO corrections to this formula are not known, but a part 



7 See, however, Ref.|78 | for a recent phenomenological application of Eq. {42}. 

IQ 



N(T), 

o.e 


ic N GBW Y=0,20,30,40 


ic N QBW Y=0,40,60,80 


; Fixed, KO . /y^ 


: ki Jr 




~_ Running, K1 / f jr 


'. ... K2 /if 


0.6 


- n « =0A iy 


- -■■ K3 if 


O.J 






0.2 


■ , 




N(x), 

o.e 


ic N° BW Y=0,20,30,40 


Y=0,40,60,80 


; K1,a =0.4 i A* 


KUcN ™ i/f 




I ... K1,a =0.2 j/ 


K1,icN MV iff 


0.6 
O.J 
0.2 




'- ~o=°- 4 if 



10"' 1 1 10"' 1 T 10 



Fig. 13: Scaling solution of BK equations ( f44b for different rapidities and initial conditions for the evolution. The 
scaling variable is defined as r = rQ sat (Y). Figure fromRef. [83|. 

of them have recently been considered in Refs. Il44l|45ll46l . In (1441 the first three terms correspond to 
the linear, BFKL, equation. They lead to an exponentially increasing scattering probability 

^linear - X~ X with A = ^^4 hi 2 . (45) 

TT 

So, the scattering probability becomes larger than 1 for small enough x, violating the unitarity bound 
which ensures probability conservation. This deficiency of the BFKL equation is rectified by the non- 
linear contribution which ensures N(r, Y) < 1 for all values of r and Y. The restoration of unitarity is 
at the heart of the saturation approaches and calls for non-linear terms in the evolution equations. These 
non-linearities correspond to the gluon fusions in the intuitive picture we provided in Section 14.2.11 
Numerical solutions of the BK equations ll83l showing the property of geometric scaling are plotted 
in Fig. [13] Similar computations including the running of the coupling Il44l |45l |46l can be found in 
Ref. Roll . A recent fit of HERA data using the BK equation with running coupling made available the 
corresponding scattering probability for computation of other observables ||84l . 

6.3 The charmonium suppression 

The first of the hard probes experimentally accessible was the J/ty suppression, proposed in 1986 by 
Matsui and Satz ||54l and first measured one year later at the CERN SPS in oxygen-uranium collisions 
Il85l . The original formulation of the problem considered the modifications of the potential (0]) in the 
case where the heavy quarks are inside a medium. The small-distance part of the potential can be com- 
puted perturbatively and the long-range Coulomb —a e g/r is screened to a short-range Yukawa potential, 
— a e ff exp{— r /r£)(T)}/r. The modification of the confining linear term in (@]) is not that simple 8 , but 
in first approximation, a deconfined medium should produce a — > 0. So, a possible modified potential 
reads (see Fig. O 

V(r) ~ -— exp{-r/r D (T)} . (46) 
r 

With this set-up, Matsui and Satz found the probability of survival for a J/$ after a thermal deconfined 
state is formed to be very small. The conclusion was that the suppression of the J/ty (as well as other 
charmonium and bottonium states) signals unequivocally the formation of a QGR 

8 This is a matter of current debate, see, for example, Ref. [ 16] and references there. 
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Fig. 14: (Left) Pictorial representation of the potential in the vacuum, Eq. ©, and the corresponding screened 
potential in the medium; (Centre) A confined cc pair in the vacuum forming a J/^P; (Right) In the case where the 
cc pair is in a thermal medium, the J is dissolved. 



The interpretation of the observed J/ ^ suppression is, however, difficult because the theoretical 
control over the modification to both hot and cold nuclear matter is not precise enough. Indeed, J/^ is 
already suppressed in hadron-nucleus collisions where the formation of a deconfined state is not possible 
owing to the small energy densities created. The charmonium states are fragile enough to be destroyed 
by interaction with the surrounding nuclear matter once produced inside the nuclei. The proton-nucleus 
programme at the SPS was essential to fix the benchmark for this normal nuclear absorption. Although 
the understanding of the J/ty suppression is still not satisfactory, a good description of the pA data is 
provided by a simple probabilistic approach in which a J/f produced at a given longitudinal position 
zq inside the nucleus travels essentially in a straight line with subsequent scatterings at z\ < zi < ... 
If the absorption or breakup cross-section for one elementary J/^— nucleon cross-section is a abs the 
corresponding survival probability at a given impact parameter can be computed in the Glauber model 



Let us take the case of pA collisions for simplicity. The J/ \P is produced at an impact parameter 
b and longitudinal position z inside the nuclei with a proton-nucleon J/ty production cross-section a^. 
Were the J not to interact with matter (as is the case for Drell-Yan dileptons, for example) the cross- 
section would be 



T 



dz Ap A (z 1 b)a x f ) = AT A (b)a^ 



(47) 



and the collisional scaling is recovered as it should be for hard processes. The J/^, however, reinteracts 
with the nuclear matter while travelling through it with an absorption cross-section, cr a bs- Repeating 
the arguments of Section |4~T1 the probability that n J/^f + N interactions take place after the J/ty is 
produced at position z is 

-i A-n 



1 



dzp A (z, 6)cr abs 



(48) 



dz p A (z,b)a ahs 

So, the survival probability corresponds to the case n = and the cross-section (|47"T ) is modified to 9 

/oo ( roc ~\ 

dzAp A (z,b)a^explj dz Ap A {z,b)a &hB \ , (49) 

where we have used again the approximation of large- A The integral in z can be done by expanding the 
exponent and using the relation 

-i n 



</:./•»(.:../)} / <h. 2 /H;- 2 .h)... j dz n p(z n ,b) = — 



dzp(z, b) 



(50) 



9 We assume here that is small to neglect terms 0(a%), i.e. only one J/^t is produced in the nucleus 
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intrinsic p with : 




Fig. 15: (Left) J suppression in pA and AB collisions for different systems and energies measured at CERN 
SPS [ 88 89 90 1 ; (Right) J/ \fr suppression in dAu collisions at RHIC as a function of rapidity compared to different 
models of J/^f production using EKS98 nuclear PDFs, figure from Ref. [91 1. 



back and forth to obtain 

= [1 - exp {-a ahs AT A (b)}} . (51) 

Cabs 

The corresponding expression for AB collisions can be obtained using similar methods to obtain, again 
in the optical approximation, 

a AB ^(b) = ^- fds [1 - e W {-a ahs AT A (s)}} [1 - exp {-a ahs BT B (h - s)}] . (52) 

^abs J 

Equations (f5TT > and (1521 provide a good description of the suppression observed at SPS in pA, 
SU and peripheral Inln and PbPb collisions. This led to the interpretation of the anomalous J/^S sup- 
pression — that which is not described by cold nuclear matter effects (l52l — as due to the formation 
of a QGP [87], see Fig. [15] At RHIC, however, where a larger suppression was expected due to the 
larger temperatures and densities of the medium, data point to a similar magnitude of the effect I92T1 . 
Other effects like shadowing have also been studied, see Fig. [T5J but the interpretation of the data is 
still under debate. One of the main difficulties is the lack of a good benchmark, both for proton-proton 
collisions, where the actual mechanism of J/ty production is not known, and proton-nucleus collisions 
where the relative magnitude of the different effects (shadowing, nuclear absorption, etc.) cannot be 
fixed by present experimental data 10 . New data from RHIC and the LHC is expected to further clarify 
the situation, including the study of other charmonium (x c , and bottonium (T, T'...) states as well 
as open charm and bottom, see, for example, Ref. [94] and references therein. 

6.4 Jet quenching 

The original idea to characterize the properties of the medium produced in hadronic collisions by the 
suppression of particles at large transverse momentum is due to Bjorken |[95l . In the simplest view, a 

10 In fact, even the validity of Eqs. l |51| > and l |52| > needs to be assessed in view of the simplicity of the probabilistic approach 
used in their derivation. In this respect, notice that the probabilistic approach can be recovered by the low-energy (incoherent) 
limit of a simplified field theoretical calculation including formation time effects |93 1. 
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partem produced in a hard process inside the medium loses energy when travelling though the produced 
matter; by elastic scatterings in Bjorken's original proposal. This energy loss leads to a strong suppres- 
sion of the high-pT particle yields due to the steeply-falling spectrum in this region of phase space; a 
10% energy loss translates into a factor of ~ 2 suppression of the yields. 

The cross-section to produce a hadron h at a given pr can be computed using Eq. (1371 ) which 
reads explicitly, see, for example, Ref. |[96l for details 

da AB ^ h s-^ f dz 2 „ , ^ 9n da l ^ kl 



dp\dy . -t^ - 

1 i,3,k=q,q,g 



/dz d& l ^~~* 
dy2—xif i / A (x 1 ,Q 2 )x 2 f j /B(x2,Q 2 ) — -~ — D k ^ h (z,^ F ) . (53) 
z dt 



The fragmentation functions D^^z, n F ) describe the hadronization of a parton k into a hadron h with 
fraction of momentum z. In the vacuum, these quantities evolve in fip with DGLAP evolution equations 
in a fashion similar to the PDFs. They are also obtained by global fits like the ones described above. 
These quantities are, however, modified in the case where this fast parton traverses a medium. A complete 
theoretical understanding of this modification is not yet known and a simple picture including only the 
energy loss of the leading particle is usually assumed; an excellent recent review on the fragmentation 
functions in the vacuum and medium is available in Ref. [971: 



D™ d (x,t) = J T ^-P( e ) J D vac (jzr e ,t) • (54) 

Here, P(e) is the probability that the particle loses a fraction e of its energy. For the moment we have 
not specified the nature of the energy loss mechanism, we have just assumed that the effect is additive to 
the usual DGLAP evolution in a probabilistic way: energy loss does not change the evolution. 

Asymptotically, the main source of energy loss of a particle traversing a medium is radiation 
l98l . In the next section we shall present the calculation of the medium-induced gluon radiation, the 
generalization of bremsstrahlung for gluons produced in a finite size medium. A main feature of the cor- 
responding spectrum is the presence of a natural collinear cut-off given by the generalization to QCD of 
the Landau-Pomeranchuk-Migdal effect. The typical singularities, which in the vacuum are resummed 
in the evolution equations, are then absent, justifying a probabilistic implementation. 

We shall not discuss here the relative importance of different energy loss mechanisms or imple- 
mentations, but just present some results taking radiative energy loss as the main source of the effect. The 
usual implementation of the effect assumes an independent gluon emission approximation, which in the 
eikonal case leads to a Poisson distribution for the energy loss probabilities — also known as quenching 
weights [991 11 



P(AE) =PoY J ^[ HdLJi dk 

k=0 ' ^ U=l ^° 



dl mcd {ui) 
diodk 




(55) 



where dl med /dujdk± is the spectrum of medium-induced gluons, the normalization is given by the prob 
ability of no-energy loss 

roc puj ^jmed 



po = exp 



du I dk^ 

duoak±_ 



(56) 



/o J o 

and e = AE/E. The total distribution contains a discrete part, the term with k = in (|55T ). giving the 
probability that the particle exits the medium unaffected, and a continuous part for a finite amount of 
energy loss AE 

P(AE) = Po 5(AE) + p{AE) . (57) 



1 A code with the numerical values is available at URL 1 101 1 . 
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All the information about the medium is contained in the medium-induced gluon radiation spectrum 
dl med /diodk± which will be presented in the next section. In the approximation in which an arbitrary 
number of scatterings with the medium are allowed, i.e., when the opacity of the medium is large there is 
a natural parameter, called the transport coefficient, q, encoding all the information of the medium except 
its total length L. These two quantities are the only two parameters of the spectrum and the goal of the 
jet quenching studies is to obtain the value of q which better describes the data once the geometry, L, of 
the medium is known. 

The question of the geometry is far from being trivial for the realistic heavy-ion collisions which 
are made in the laboratory. Although there is a rather good experimental control on the centrality of the 
collision, the actual energy density profiles of the created medium still present some uncertainties. For the 
jet quenching phenomenology, the implementation of a geometry as realistic as possible translates into 
a better determination of the transport coefficient q, or the corresponding parameter in other approaches. 
An example of this is the determination of this parameter performed within the same approach but using 
different geometries: 

1. Taking a fixed length of L = 6 fm (roughly the average distance traversed by a parton produced 
inside a gold nucleus of radius L ~ 1.2A 1 / 3 and travelling perpendicular to the beam direction) 
the estimated value is q ~ 1 GeV 2 /fm [100]. 

2. Taking the same geometry, a cylinder of radius L = 6 fm, but now computing the suppression of 
the yields before averaging in the length, the resulting value is 12 q ~ 5 4- 15 GeV 2 /fm 11021 . 

3. Taking a geometrical profile from a hydrodynamical description of the heavy-ion collision the 
obtained value is q ~ 8 GeV 2 /fm 111041 . Now, the medium is not static but evolving in time with 
longitudinal and transverse expansion. In this case, the average value of the transport coefficient 
depends on the trajectory — each different path leads to a different average of q. 

In recent years an attempt was started to fix the medium properties by common fits to several sets 
of data lfT^[Totl[T07l[TMl[TMl[T04l . 

6.4.1 Hydro meets jet quenching 

The hydrodynamical description of the soft part of the spectrum in heavy-ion collisions makes this ap- 
proach very appealing as a realistic geometry implementation of the energy density profiles for the jet 
quenching calculations. The usual approach is to define a local relation between the transport coefficient 
and one of the properties of the medium which can be accessible by a hydrodynamical model like the 
energy density lfTT0lfT071[TM[T09l[T04l 

q(x,y,T) = 2Ke 3 /\x,y,T). (58) 

Here, for simplicity, it has been assumed that the particle direction is perpendicular to the beam axis. 
The factor of 2 in (|58T > is a useful convention, as, taking K = 1, the estimate of the transport coefficient 
for an ideal QGP is recovered H 1 X OB 1 3 . 

The results from Refs. H X 1 1 1 11001 11121 1 1 1 311 allow us to compute for a dynamical medium (in 
which the transport coefficient decreases as a function of time due to the longitudinal and transverse 
expansion of the medium) an equivalent static scenario in which the transport coefficient is given by 

# (£ " ft) <l(0 , (59) 

Using this relation and the local values of the transport coefficient defined by d58l ), the suppression of 
high-pT hadrons can be computed using the factorized expression for the cross-section (|53T > and the 

12 The same value is obtained if instead of a cylinder a Woods-Saxon distribution is assumed for the energy density 1 103 1. 
13 Notice that although this is only an estimate, all perturbative estimates lead to similar numerical results for similar medium 
temperatures, a s , etc. so, we take the if-factor as a departure from a perturbative estimate of q at leading order. 
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Fig. 16: Nuclear modification factors Raa for single-inclusive and Iaa for hadron-triggered fragmentation func- 
tions for different values of K — 2, 3, 4... The best fit is marked by a wider green line. 



corresponding medium-modified fragmentation functions (154] ) calculated using the quenching weights 
(1551 . In this set-up, once e(x, y, r) is known, the only free parameter is then K to be fitted to experimental 
data. The following results take these energy density profile functions from the hydrodynamical results 
in Ref. 11141 : a code with the numerical output is publicly available at URL HI 15 II . 

The first data to be taken into account is the suppression of light partons; we take into account only 
7T° in the analysis to avoid contributions from baryons. The suppression is normally plotted in terms of 
the nuclear suppression factor 

da AA /dp 2 T dy 

KaA = ITt 3 / , 2 , • ( 6 °) 

In Fig. [16] the results for the suppression of 7r°'s produced at high p? in central AuAu collisions 
measured by the PHENIX Collaboration at RHIC II 1061 are plotted together with the theoretical calcu- 
lations for different values of K. Highlighted is the value corresponding to K = 4 which is the one 
providing the best description for this particular set of data. 

In order to obtain additional constrains to the value of K (or q), other observables can be studied. 
Increasing interest at RHIC is put on the suppression of back-to-back correlations and also on heavy 
quarks. The latter will be commented on later. The two-particle correlations, on the other hand, are 
expected to be more sensitive to the medium properties as the surface bias effects would be smaller than 
in the one-particle inclusive case. In a two-to-two perturbative process like the ones considered here, 
the two initial partons are produced back-to-back and hadronize independently. This corresponds to a 
different modification of the fragmentation functions for each of the partons, as they follow different 
path lengths through the medium. There is, however, no extra free parameter in the calculation. In 
the insertion of Fig. [16] the suppression of the back-to-back azimuthal correlations measured in central 
AuAu collisions by the STAR Collaboration at RHIC [116] are compared with the theoretical curves for 
the same values of K. This is a check that these two sets of experimental data can be described with the 
same medium parameter K ~ 4. The corresponding definitions are 

D aa {z t ,Pt 1S ) ~ n 

Dpp{ZT,P T ) 



where the hadron-triggered fragmentation function is defined as 111 171 



D aa (z t ,Pt 18 ) = Pt 



da^ 2 /dy tr ^dp% lg dy &ssoc dp^ 
da h A \/dy^dp^ 



.assoc 



(62) 



Here, = p^ soc / p^ lgs . It should be noticed that the objects defined in this manner are, in general, 
very different from the fragmentation functions measured, for example, in e + e~ annihilations at similar 
virtualities owing to the strong bias produced by triggering in a steeply falling perturbative slope. 

The average value of the transport coefficient obtained is, then, (q) ~ K2(e) ~ 8 GeV 2 /fm. 
6.4.2 The heavy quarks 

Gluon radiation is suppressed by mass terms in the vacuum by a mechanism called the dead-cone effect. 
These mass terms appear in the propagators of the quarks before the radiation point and the corresponding 
suppression factors are given by the ratio of the massive over massless propagators E118I1 



Here, the dead-cone angle 9\ c ~ M 2 / E 2 defines the angle below which radiation is suppressed. This 
suppression can be understood intuitively by the inability of the heavy quark to build coherently the wave 
function of the gluon due to the different velocities of massless and massive particles; in the limiting case 
6 — > the radiation should cancel as no constructive interference is possible. 

In the case of the medium 111 19111201 11211 [1221. the radiation is also suppressed by terms similar 
to (l63l in the propagators, but, interestingly, the LPM suppression is now less effective as the gluon 
formation time is smaller for a massive quark than for a massless quark: the suppression of the radiation 
at small angle is smaller, so the radiation is larger. The net effect is a competition between the suppression 
due to mass terms in the quark propagators and the smaller suppression of collinear gluons due to LPM 
interference. However, although in some limited regions of phase space the dead cone is filled II122L 
the net effect is a suppression of the radiation by mass terms. So, the prediction from the formalism 
is a smaller energy loss for heavy quarks than for light quarks. This different effect is controlled by 
mass terms in the medium-induced gluon radiation expressions and implies no new free parameter in 
the calculations, providing additional constrains to the determination of the medium properties and the 
nature of the energy loss. 

Recent attempts to measure the jet quenching of heavy quarks have been made at RHIC 11231 
11241 11251. Without a proper vertex reconstruction, the heavy quarks are measured by the semileptonic 
decays of heavy hadrons into electrons after subtraction of those coming from photons. These non- 
photonic electrons originate from the decay of both charm and beauty hadrons and the experimental 
access to the relative contributions is rather indirect. The theoretical situation is no better and a large 
indetermination of this quantity is present. State-of-the-art calculations of heavy-quark cross-sections 
include resummations of large logarithms originated by the mass of the heavy quark in a scheme called 
Fixed-Order-Next-to-Leading-Log (FONLL) approximation 11126111271 . In Fig. [F7]the comparison of the 
non-photonic electrons coming from heavy hadron decays into electrons H 12811 with experimental data in 
proton-proton collisions at yfs = 200 GeV at RHIC is shown. The agreement is reasonable within the 
theoretical uncertainty band computed by varying quark masses and renormalization scales, although the 
experimental data have the tendency to be underestimated. The relative contributions from c and b quarks 
are also plotted for the extremes of the bands. When computing the suppression due to medium effects 
[129], these uncertainties in the proton-proton benchmark translate into a corresponding uncertainty 
band for Raa as mass effects are more important for beauty than for charm. This calculations have been 
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Fig. 17: (Colour online) (a) Left: Comparison of the FONLL calculation of single inclusive electrons II 1281 to 
data from pp collisions at y/s = 200 GeV 11231 . Upper and lower lines are estimates of theoretical uncertainties, 
obtained by varying scales and masses, for details see the text, (b) Right: The nuclear modification factor R% A 
of electrons in central Au-Au collisions for an opacity of the produced QCD matter characterized by the time- 
averaged BDMPS transport coefficient q = 14GeV 2 /fm. The shaded band indicates the theoretical uncertainty 
of the perturbative baseline only. Red dashed and blue dotted curves show R e AA for c-quark and 6-quark decay 
contributions, respectively. Data taken from Refs. 11124111251 . Figure from Ref. 11291 . 

performed with a Woods-Saxon profile for the medium, i.e., a static case, but the main conclusions hold 
equally for the case of a hydrodynamical medium profile. 

This exercise shows the need for a more precise identification of the charm and the beauty mesons 
for a correct interpretation of the dynamics underlying energy loss of heavy quarks. It is worth mention- 
ing that the mass effects in the electrons coming from decays of charm hadrons alone are rather small 
and the suppression is almost mass-insensitive for > 5 GeV. The effect is much more important for 
beauty. So, although taken at face value, the suppression in the experimental data looks to be underesti- 
mated by this analysis, a larger contribution of the charm quark would help in reducing this discrepancy 
1122111201. 

It is also worth mentioning that the formalism explained here is likely to be pushed too far, spe- 
cially for beauty, and other effects could appear. One of the extra effects which have been advocated in 
recent years is a large contribution of the collisional energy loss which would not be negligible for RHIC 
kinematics IU311ll32lll331ll34ll . Other possibilities are a larger contribution of the heavy baryons (which 
have a smaller branching ratio to electrons, explaining part of the suppression) in nuclear collisions at 
intermediate transverse momentum H135L 11361 : or non-perturbative modifications of the hadronization 
of heavy quarks, which owing to the smaller hadronization times and Lorentz factors would hadronize 
inside the medium 1113711 . 

With present experimental data, none of the above descriptions is completely satisfactory and a 
more precise determination of the mass effect in the heavy-quark suppression would be needed. This 
identification will be possible with the detectors of the different LHC experiments and also with the 



71 



upgrades planned for the RHIC experiments. 



6.5 New developments in jet studies in heavy-ion collisions 

A jet is a bunch of particles flying to the detectors in a given limited region of the whole phase space 
available. This experimental definition corresponds to the production of high-p-r quarks or gluons in hard 
interactions or in the decay of a very massive particle. In general we can associate large virtualities to 
these produced partons which then radiate (mainly gluons) to become on-shell and eventually hadronize. 
So, jets are extended objects corresponding to the large difference in the scales involved: the large pr, 
defining the typical virtuality of the object, and the hadronization scale O(Aqcd)- hi other words, 
the transverse momentum of the produced parton determines the maximum phase space available for 
radiation. 

The characterization of jet properties in the vacuum is one of the best tests of QCD. The most 
favourable experimental condition correspond to e + e~ annihilation into hadrons. Jets in lepton-proton 
or proton-proton collisions need a better control over the background. This is specially relevant for the 
LHC where the background for the underlying event and/or multiple collisions per crossing in proton- 
proton collisions will be important. The background will be even larger for central PbPb collisions; 
peripheral or minimum-bias PbPb collisions will be similar to proton-proton collisions in this respect 
due to the reduced luminosity for the formed and the slightly smaller centre-of-mass energy. 



6.5.1 The vacuum 

The description of jets in the vacuum is a very well developed subject with a large amount of literature, 
see, for example, Ref. If 1381 . Here we only provide some heuristic arguments based on the LO descrip- 
tion of the effect which is specially relevant for the implementation in Monte Carlo event generators such 
as PYTHIA lfl39l . HERWIG lfl4tjl . or SHERPA lfl4TI 

The spectrum of gluons radiated with a given energy uj and transverse momentum kr by a highly 
energetic quark or gluon of energy E produced in a hard collision is 

dl™ = «« J_ p( x ^ J_ 1 m 

dzdk 2 ^ 2tt k\ z ~* 2ir kj, z ' 

where kx is the transverse momentum of the gluon with respect to the original parton direction and 
z = uj/E the fraction of momentum carried by it after the radiation. Equation (|64l ) contains infrared 
(z — > 0) and collinear (9 ~ k^jui — ► 0) divergencies. The available phase-space to integrate the 
radiation is determined by the virtuality Q 2 ~ p\ as mentioned above, introducing a cut-off, [j? for the 
collinear divergency: 



jrvac 

dk T-^T ~a s log(Q7^). (65) 



So, even when a s is small, the presence of large logarithms from extended integration regions in phase 
space makes these higher order terms important and a resummation of an arbitrary number of radiated 
gluons is necessary. This leads to a gluon emission which is ordered in virtuality. 

An intuitive way of understanding the formulation of the problem is through the radiative decay 
analogy HI 391 in which: i) there is an ordering variable t — for a radiative decay it is the time, in our 
case it is the virtuality; ii) the probability that something will happen at "time" t (a decay or a parton 
branching) is given by a function f(t); iii) there is an additional requirement that something can only 
happen at time t if it did not happen at earlier times. In these conditions, the probability that nothing has 
happened, M(t), is given by 

^ = -f{t)N(t) => N(t) = exp j- jf f(t')dt'} . (66) 
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So that, the probability that there is a branching or a decay at time t is just 



V(t) = f(t)exp^- 1^ f(t')dt'Y (67) 

In our case of a partem branching, f(t) is known and the resulting probability to radiate one gluon at a 
given virtuality Q 2 from a parton produced with virtuality Q 2 nax would be (see also Ref. [138] for more 
details) 

dv(z, q 2 ) di™ c f r Q ™°* , r , d/ vac 1 d/ vac . ,„ 2 N . 

^ i =^ exp {-y Q2 ^^l^^w^^- (68) 

Notice that removing the virtuality ordering in (I68T ) would lead to a Poisson distribution for the multiple 
gluon radiation of the form (I55T ). Writing the spectrum explicitly, 



dV(z,Q 2 ) = ^L-Lp( z )A (QLx,Q 2 ) dzdQ 2 . (69) 



The term A(i), defined as 



A(i) = exp 



J to * J ^ 



(70) 



is called the Sudakov form factor and gives the probability that no radiation took place between the two 
scales £ Q and i. This approach can also be applied to the initial-state radiation for the PDFs. Now the 
virtuality ordering goes in the opposite direction, increasing from an initial Qq to a given Q 2 value. So, 
if the initial condition for a given PDF is f(x, Qq), its value at Q 2 is given by the iteration of (l69l 



f(x,Q 2 ) = A(Q 2 )f(x,Q 2 ) + f f *Zp(z)f (- z ,t) . (71) 

The first term is the contribution in which no resolving radiation corrects the initial value of the PDF 
between the two scales and the second one gives the evolution due to radiation. Equation (TTTb is the inte- 
gral representation of the DGLAP equations for only one flavour, which in its more familiar differential 
form stands as 14 

df(x,Q 2 ) a s f dz fx 



d\ogQ 2 2tt 



z 



P(z)fl-,t). (72) 



So, the DGLAP evolution equations take care of the logarithms d65l ) and the same evolution applies for 
PDFs and jets in the appropriate kinematical regimes. 

The probabilistic approach given by Eq. (TTTb is specially interesting for Monte Carlo implemen- 
tation. The parton shower algorithms in all of these codes follow the ideas sketched above with some 
modifications for each particular implementation. We can, in fact, understand these processes as an al- 
most independent gluon emission approximation, in which the presence of an ordering variable breaks 
the complete independence. The ordering variable we have been presenting above is the virtuality, but 
other more or less sophisticated ordering variables are also used. For example, HERWIG defines a vari- 
able which leads to an ordering in the emission angle of the subsequent partons. This angular ordering 
encodes colour coherence effects which are absent otherwise. 

The Monte Carlo implementation of the parton shower in event generators follows an iterative 
procedure; for definiteness, here we take the case of a time-like (final-state) evolution of a produced jet. 
The basic problem consists in generating the virtuality t-i and fraction of momentum X2 for a branching 
which takes place after a branching with virtuality t\ and fraction of momentum x\. 



14 The splitting functions in the integral form are, however, regularized through the + prescription whose origin is the virtual 
corrections when computing the one-gluon emission, see, for example, Ref. 11381 . 



Fig. 18: Result of the parton shower algorithm explained in the text for the hadronization of a quark-antiquarkpair 
produced in the decay of a virtual photon. The colour flow in only one of the hemispheres is depicted as colour 
lines. 



1. First, the virtuality ^ is decided by choosing a random number 1Z and solving 

A(t2) 



A(ti 



n. (73) 



This equation makes use of the probabilistic interpretation of the Sudakov form factor explained 
above. 

The next step is to compute the value of the fraction of momentum of the radiated gluon (22) by 
solving 

fXi/xi rl-e 

dz^-P{z)=U ! J dz^P(z), (74) 

with VJ another random number and e an infrared cut-off. 

3. Once the values of ti and X2 are obtained, checks are made; these could depend on the actual MC 
implementation. If the branching is possible, a new branching is tried by repeating points 1, and 
2, above with the appropriate kinematical variables. The branching stops when some cut-off scale 
to = 0(1 GeV) is reached. 

4. The produced partons (quarks and gluons) are colour-connected in a definite manner so that colour- 
less objects can be found. These colourless objects are hadronized by some non-perturbative model 
into final particles. The hadronization model depends on the actual MC implementation, for exam- 
ple PYTHIA assumes that strings are formed and decay by subsequent breaking, while HERWIG 
forms colourless clusters of quarks and antiquaries which then decay into the final particles. 

The typical structure of this shower iterative procedure can be seen in Fig. [18] where a virtual 
photon produces a qq pair in opposite directions, each of them forming a parton shower. The colour 
connection between the different partons is depicted as colour lines for only one of the hemispheres. 

The description of experimental data provided by the main event generators is excellent in those 
quantities sensitive to the parton shower structure. For this reason, these codes are essential tools in 
the phenomenological interpretation of the data and constitute, in fact, one of the main parts of the 
experimental software for data analysis. 
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6.5.2 The implementation of the medium effects 

The vacuum radiation process is a very mature subject in which a large amount of work both from 
theory and experiment has led to a very good knowledge of the effects involved. The counterpart for the 
medium radiation is, in comparison, in its infancy. Indeed, on both the theoretical and the experimental 
side, this is a much more recent subject where a lot of progress is expected in coming years. A learning 
period with dialog between theory and experiment will be needed. A breakthrough in this respect is the 
first measurement of reconstructed jets in heavy-ion collisions by the STAR Collaboration 111421 which, 
although suffering from severe experimental uncertainties, paves the road to future progress in the field. 

Several approaches to the medium implementation of parton showers exist at present H1431 11441 
ll451U46l [T47i. none of them being completely satisfactory but addressing different issues in the problem. 

There are several effects which could modify the evolution of a parton shower in a medium after 
the hard process, in particular: 

1. Medium-induced gluon radiation modifying the vacuum radiation structure. As we have seen, this 
radiation is free of IR or UV divergencies, but is enhanced by powers of the traversed medium 
length a s L. The angular structure is also modified. 

2. Non-eikonal corrections to the parton propagation in medium (also known as collisional energy 
loss) 

3. Energy flow from and to the medium — which could eventually modify the hydrodynamical evo- 
lution of the medium. 

4. Modification of the non-perturbative hadronization. If the transverse momentum of the particles is 
large enough, the non-perturbative (slow) processes would take place in the vacuum, so, as a first 
approximation all models assume vacuum hadronization. 

5. Modification of the colour structure of the shower evolution. I.e., the colour structure plotted in 
Fig. [18] corresponds to the vacuum, where colour exchanges are only possible within the created 
particles. In the case where this structure is formed in the medium, colour exchanges with the 
particles in the medium are also possible, leading to a modified colour structure of the parton 
shower. 



One of the main issues to be solved is to find the appropriate ordering variable in the presence 
of a medium. For the vacuum, the space-time picture of the parton shower is, in fact, irrelevant as no 
external scale exists. The spatial extension of the medium indicates, however, that time should play a role 
as an ordering variable; there are, in fact, known interferences between the length of the medium and the 
radiation, leading to LPM suppression, which are included in the computations of the medium-induced 
gluon radiation. 

A simple implementation of medium effects has recently been proposed which attempts to include 
in a consistent manner the vacuum and medium contributions to the parton shower 111461 1 1471 1 1481 [T49T1 . 
In this approach, a medium term is included in the splitting function 



-tot 



(z) = P vac (z) P tot (z) = P vac (z) + AP(z, t, q, L, E) , (75) 

as given by the medium-induced gluon radiation spectrum by matching the vacuum term, recalling the 
definition of the medium-induced gluon radiation as the total radiation in the presence of a medium with 
the vacuum contribution subtracted: 

dI u d/VaC + u d/m6d (76) 
du dk± du> dk± duo dkj_ 

The vacuum spectrum defines the (vacuum) splitting function, which at small fraction x = 1 — z is 

tlP '' H ^^P vac (z), P™(z)^^-. (77) 



dzdtf 27Tk2 ^ h ^ ' - 1 
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Fig. 19: Sudakov form factors for quarks (upper curves and symbols in each plot) and gluons (lower curves and 
symbols in each plot), for parton energies E — 10 GeV (upper plot) and E = 100 GeV (lower plot), and for 
different medium parameters. 



This indicates the following matching for the definition of the medium-modification term of the splitting 
function lfl48l [1491 

, . 2vrt dl med 

AP(z,t)~ — — , (78) 

a s azat 

where u = (1 — z)E and = z(l — z)t are taken in the medium-induced gluon radiation spectrum; 
parton masses have been neglected. 

The new splitting functions are also implemented in the calculation of the Sudakov form factors, 
the main component in the parton shower algorithms. The corresponding modifications are plotted in 
Fig. [I9]where the vacuum and medium cases are compared for different values of the medium properties 
encoded in the jet quenching parameter q. This comparison clearly shows the enhancement of the radi- 
ation expected by this new term in the splitting probability which will be translated into larger intrajet 
multiplicities and different jet structures. 

In general we expect the jets produced by the modified parton shower to present: 

1. Larger multiplicities due to the enhancement in the radiation probability 

2. Softer spectra due to the energy loss of the partons by radiation, in particular that of the leading 
parton 

3. Larger radiation angles (jet broadening) due to the absence of collinear divergencies in the medium 
radiation contribution. 

These general expectations are independent of the particular implementation of the effect and are solely 

1? 



based on the properties of the medium-induced gluon radiation spectrum. A Monte Carlo implementation 
allows one to study the effect of important kinematic constrains which are otherwise difficult to address 
in analytical calculations. 

In order to take into account the interplay between the medium length and the virtuality (or angu- 
lar) ordering, a mixed approach is used in Ref. 11461 [1471. For each splitting, the formation time of the 
gluon is computed as 

^form — "To" (79) 

and the length in the splitting function (l78l) for the next branching corrected for this quantity. This 
effectively imposes a veto for radiation during the time in which a branching is being formed. This time 
is taken from the average value d79l ). So, although the ordering variable is the same as in the vacuum, 
some medium-radiation is vetoed by these formation time arguments. 

The modified splitting probability and the corresponding corrections for length evolution are im- 
plemented as modifications of the shower routines in PYTHIA and HERWIG and available for public 
use Ill461ll47lll50[|15fi . These are not official releases of the corresponding codes but just modifications 
of them. The nicknames Q-PYTHIA and Q-HERWIG are used for the corresponding implementations 
11461 IMZl [BO). 

The intrajet parton distributions computed with Q-PYTHIA are plotted at the parton level (before 
hadronization) in Fig. |20]for the fragmentation function (as a function of £ = — log(z)); the transverse 
momentum with respect to the jet axis; and the angle with respect to the jet axis. The main expectations 
listed above, the softening and enhancement of the multiplicity and the jet broadening, are clearly visible 
in this plot. 




Fig. 20: Intrajet parton distributions in £ = log(£j e t/p) (left), px (middle), and 8 (right) for a gluon of initial 
energy Ej et — 100 GeV in a medium of length L = 2 fm and for different transport coefficients q = (black), 5 
(red), and 50 (blue lines) GeV 2 /fm. 



7 The parton propagation in matter and the medium-induced gluon radiation 

In this section we provide general ideas on the calculation of the medium-induced gluon radiation and 
some of the properties of the spectrum. This spectrum has been computed by several groups 111521 11531 
I154[|155l fl561 using different formalisms and approximations. For a more complete derivation, we refer 
the reader to the available lectures Il60ll and reviews |[57ll58ll59ll , or the original papers. 

To describe the jet quenching phenomenon, we start from a high-px quark or gluon produced 
in an elementary hard collision which subsequently interacts with the surrounding matter. So, the first 
question we need to address is how a highly energetic particle propagates through a dense medium. 
In the case of very energetic particles, a semiclassical approach is possible in which changes in the 
medium configuration due to the passage of the fast particle, recoil, are neglected. In these conditions, 
the propagation of a particle at transverse position x to transverse position y from time t\ to can be 
described by the Feynman path integrals with a background field A(x) (see Fig. l2Tb : 



G(x,ti;y,t 2 ) = / £>r(t)exp 



E 



d£ 



+ ig / d£A(r(0) 



(80) 
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Fig. 21: Two paths as described by Eq. (f80b . Figure from Ref. 11571 . 

In the asymptotic case, E — » oo, only the paths in which dr/d£ = survive, i.e., the particle does 
not change its transverse position, but travels in a straight line. In this case, the propagation is given by 
the Wilson line 



W(x) = exp Ug J d£A(x) > . (81) 

Now, the propagation of a quark or a gluon through QCD matter (described by the fields A(x)) 
changes the phase of its wave function by rotation in colour space. 

To continue, we notice that any observable will involve colourless combinations of the Wilson 
lines, which implies that, at least, the colour trace of two Wilson lines should be preset for the colour to 
cancel. This colourless object will, then, measure the state of the medium in a given configuration. To 
compute a given observable, an average on the ensemble of possible medium configurations needs to be 
performed. So, the simplest object which will appear in any calculation is 

' Tr(W(x)W(y)) , (82) 



iV 2 -l 



where the prefactor l/(iV 2 — 1), with N the number of colours, has been included to average over the 
initial colour configurations — in this case for a gluon, for a quark the prefactor would be \/N. This 
object contains all the non-perturbative physics and, in the calculations presented here, only combinations 
such as (l82l appear. 




Fig. 22: The medium-induced gluon radiation diagram 



One of the main issues in this type of analysis is to compute the medium averages (1821 ) for which 
several prescriptions exist. For a medium in which a large number of scattering centers interact with the 
propagating parton, an opaque medium, a widely used approximation is 



J J- l Tr(W( X )W(y)) ~exp{-~ J d^(0(x-y) 



(83) 



This defines the transport coefficient, q, as the pref actor of the typical small-distance r 2 -dependence 
of QCD dipole cross-sections. This prescription corresponds to propagating partons which describe 
Brownian motion in the transverse plane characterized by q. The transport coefficient can, hence, be 
interpreted as the average transverse momentum squared per mean free path: 



A 



(84) 



For the jet quenching calculations in the previous sections, we need to compute the medium- 
induced gluon radiation. This implies the propagation particle with large energy E, let us say a quark, 
which radiates a soft gluon with energy uj <C E at a small angle 6 ~ k±/u. The typical diagram to be 
computed can be seen in Fig. |22] where the blobs represent scattering with the background field. Figure 
|22]contains the propagation of three particles in the medium, the initial quark and the produced gluon and 
quark. Each of these propagations are described by (l80l ) — in fact, the quark is considered completely 
eikonal, E — ► oo, so that the Wilson line (I8TI ) describes the propagation. In these conditions, the only 
non-perturbative object is the average of the Wilson line (l82l . The final result takes a very compact form 



dl 
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where 



K, (r(x),x; r(x),x\u) 



f \ v r 

J Prexp i— J d£ 



dr 



2uj 



(85) 



(86) 



corresponds to a 2-dimensional harmonic oscillator with time-dependent imaginary frequency. In Eq. 
d85l the three terms correspond to (i) the gluon emitted inside the medium in both amplitude and conju- 
gate amplitude; (ii) the emission inside the medium in amplitude and outside the medium in conjugate 
amplitude; (iii) and when the gluon is emitted outside the medium in both amplitude and conjugate 
amplitude (see Fig. l23l . 

The last contribution corresponds to the vacuum radiation which is normally subtracted to define 
the medium-induced gluon radiation as 

dl med dr &c 

+ cj , ;91 . (87) 
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Fig. 23: The three contributions to the squared amplitude of the medium-induced gluon radiation. The dashed line 
is the cut indicating the final outgoing particles. 




Fig. 24: Left: numerical results for the medium induced gluon radiation spectrum uidl mcd / 'duidk± [Eqs. (f85T > and 
(|87T >1 of a quark in a static medium as a function of the dimensionless variables (l88l l. Right: Same but integrated 

in k± < to 



The spectrum d85l > with the vacuum subtraction d87l ) is the building block of the calculations presented 
in Section \6A\ 

7.1 Heuristic discussion 

In Fig. [24] we present the results for the double-differential medium-induced gluon radiation spectrum 
for a quark traversing a static medium of length L. The results are given as a function of the variables 

ui c = \qL\ k 2 = %. (88) 
2 qL 

One important feature of the spectrum is the presence of small-fc^ and large-w cuts which can be under- 
stood by the formation time of the gluon, 



^■forrn — ,2 i (89) 



2lo 

which is controlled by the dynamical term in the path integral d80l ). To see this we can expand the path 
integral in the number of scatterings 

G(x,t ;y,t/) = G {t -► t f ) + G (t -► t^igA^Goih -►*/) + 

+^G (t -» ti)i 5 A(xi)G (ti -» t2)*^(x a )Go(t2 -»*/)... (90) 
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where Go(U — ► t,) is a short-cut notation for the vacuum propagator from time ti to tj with transverse 
positions defined by Xi and xj, respectively. The vacuum propagator can be written as the Fourier 
transform of a plane wave 

Go(x i± ,t i; x (i+1)±) i (m) ) = J |^ e ^(fe-^ + i)) e -Pix(x i x-x (i+1)x ) (91) 

In d9Tb the exponential factors define the coherence length d89l ): when tf orm 3> tj — the phase is 
small and we can write 

d 2 p. 



nf d Pi± P -ipi±( x i-L- x (i+i)±) [fa . p-«Pi±( x -y) 
. J (2n) 2 ^ ./ dX±e 



(92) 

So, all the scatterings d90l) take part at the same time — coherently. On the opposite side, when the for- 
mation time is much smaller than typical lengths in the medium, the exponentials oscillate very fast and 
only the first term survives. In other words, in the incoherent case, the cross-section counts the number 
of scattering centers, while in the coherent case the whole medium acts as a single scattering center. As 
a result, a reduction of the gluon radiation is produced in the latter case. This is the generalization to 
QCD of the Landau-Pomeranchuk-Migdal effect. The numerical effect appears clearly in Fig. |24]as a 
suppression of the spectrum for small values of k 2 . An important consequence is that the spectrum is 
neither collinearly divergent (i.e., it can be safely integrated to k± = 0) nor infrared divergent (i.e., it 
can be integrated to uj = 0) as can be seen in Fig. [24] In contrast the vacuum part of the spectrum (l85T ) 
presents both collinear and soft divergences, 

dl vac a s C R 1 1 

(93) 



dujd 2 k± it 2 k 2 ^_ uj 



The position of the infrared cut in the k± -integrated spectrum (Fig. |24] Right) can also be understood in 
terms of the formation time: integrating the spectrum in the kinematically allowed region < k t < uj, 
and noticing 15 that (k 2 ) ~ y/gu then a suppression of the spectrum for u < q 1 ^ 3 should appear. 
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